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Aims. The main goal of the present study is to determine the fractional SiO abundance in high-mass star-forming cores, and to investigate its 
dependence on the physical conditions. In this way we wish to provide constraints on the chemistry models concerning the formation of SiO in 
the gas phase or via grain mantle evaporation. The work addresses also CH 3 CCH chemistry as the kinetic temperature is determined using this 
molecule. 

Methods. We estimate the physical conditions of 15 high-mass star-forming cores and derive the fractional SiO and CH3CCH abundances in 
them by using spectral line and dust continuum observations with the SEST. 

Results. The kinetic temperatures as derived from CH3CCH range from 25 to 39 K, the average being 33 K. The average gas density in the 
cores is 4.5 • 10 6 cirr 3 . The SiO emission regions are extended and typically half of the integrated line emission comes from the velocity 
range traced out by CH3CCH emission. The upper limit of SiO abundance in this 'quiescent' gas component is ~ 10~'°. The average CH3CCH 
abundance is about 7 • 1(T 9 . It shows a shallow, positive correlation with the temperature, whereas SiO shows the opposite tendency. 
Conclusions. We suggest that the high CH3CCH abundance and its possible increase when the clouds get warmer is related to the intensified 
desorption of the chemical precursors of the molecule from grain surfaces. In contrast, the observed tendency of SiO does not support the 
idea that the evaporation of Si-containing species from the grain mantles would be important, and it contradicts with the models where neutral 
reactions with activation barriers dominate the SiO production. A possible explanation for the decrease is that warmer cores represent more 
evolved stages of core evolution with fewer high-velocity shocks and thus less efficient SiO replenishment. 
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1. Introduction 

High-mass star-forming regions with large molecular column 
densities and bright molecular lines offer excellent opportu- 
nities to study interstellar chemisty. Results from chemistry 
models can be combined with observational data on the phys- 
ical conditions and molecular abundances. This information 
can help us to link the properties of massive dense cores, 
their chemical characteristics, and phenomena like massive 
outflows, ultracompact (UC) H 11 regions and molecular masers 
into a coherent evolutionary sequence. This is useful since 
the process of high-mass star formation is not yet well un- 
derstood (see, e.g., Garay & Lizano 1999| IFontani et al. 2 002 
[Thompson and Macdoira!d2003| [ 

In this paper we derive the abundances of silicon monoxide, 
SiO, and methyl acetylene, CH3CCH, in a sample of massive 
molecular cloud cores, and discuss the relation of these abun- 
dances to the physical conditions and the probable evolution- 
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ary stages of the cores. The purpose is to contribute towards 
a better understanding of the interplay between the processes 
associated with massive star formation and the physical and 
chemical properties of the surrounding dense core. 

SiO is believed to trace exclusively shocked gas as 
its spectral lines have wings and are often shifted rel- 
ative to the emission fr om the ambient gas (see, e.g., 
iMartin-Pintado et al. 19921 ISchilke et al. 1997i . Current un- 
derstanding is that SiO is evaporated from the dust grains 
when the shock velocity is greater than about 20 km s _1 . 
This, however, depends on the grain-mantle composition in 
the preshock gas (e.g., ISchilke et al. 1997l l. In dense and 
warm giant molecular cloud (GMC) cores, SiO seems to 
present also in the quiescent gas component. In few inter- 
ferometric maps available extended SiO emission at ambient 
cloud velocity is either seen in the vicinity of high-velocity 
outflows ( Hat chell et al. 20(711 or totally separate from them 
(ILefloch et al. 19981 ICodella et al. 1999l|Shepherd et al. 2004| 
IFuente et al. 20 05 1. The fractional abundance of 'quiescent 
SiO' is estimated to be about 10" 10 in GMC cores. SiO ab- 
sorption measurements towards dormant GMCs, so called 'spi- 
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ral arm clouds', give similar abundances (Grea ves et al. 19 96 1. 
The present observational evidence is insufficient to determine 
whether the presence of SiO is caused by thermal evapora- 
tion, enhanced neutral-neutral production pathways, photon- 
induced reactions, or by shock removal. To resolve this indis- 
tinctness it seems beneficial to determine fractional SiO abun- 
dances and kinetic temperatures in a representative sample of 
GMC cores, and to correlate these also with the kinematic in- 
formation contained in the spectral lines. 

CH3CCH is an organic molecule observed rather frequently 
towards dense cores. Its abundance is on the order of 10~ 9 . 
The rotational temperature, T rot , derived from a series of 
Jk — (J - rotational lines of CH3CCH is considered a 
good estimate of the gas kinetic temperature, 7\j n , in molec- 
ular clouds (e.g.. lAskne et al. 19841 |Bergin et al. 1994| >. Since 
CH3CCH has a relatively low dipole moment (pi = 0.78 D, 
Bauer etal. 19791 iBurrell et al. 19 80 1, its rotational levels are 
thermalized at densities n(H2) ~ 10 4 cirT 3 or higher (e.g., 
|Kuiperetal. 19 84 1. This property together with the assumption 
that CH3CCH emission is optically thin (r <K 1), allow us to as- 
sume local thermodynamic equilibrium (LTE), so that Tu n may 
be derived. Besides for 7\i n , we use T mt for estimating the dust 
temperature, Tj, which is needed to derive the total gas column 
densities (~ N(H2)) from dust continuum observations. 

The source list consists of 15 high-mass star-forming cores 
associated with OH, H2O, and CH3OH masers, UC Hn re- 
gions and bright, thermal SiO rotational line emission. The 
sources were selected from the SiO survey of |Harju et a~1 998 
(hereafter HLBZ98). The positions selected for the SiO and 
CH3CCH line observations are listed in Table |2] These posi- 
tions served as map centres in the continuum observations. 

The SiO and CH3CCH column densities and the kinetic 
temperatures in the present paper are derived from spectral line 
observations with the SEST. The estimates for the molecular 
hydrogen column densities, N(&2), and the core masses are 
derived from dust continuum maps obtained with the SIMBA 
bolometer at SEST. The observations and the data reduction 
procedures are described in Sect. 2. The direct observational 
results are presented in Sect. 3. In Sect. 4 we describe the meth- 
ods used to derive the physical/chemical properties of the GMC 
cores. In Sect. 5 we discuss the results, and finally, in Sect. 6 
we summarize our major conclusions. 

2. Observations and data reduction 

2.1. Molecular lines 

The spectral line observations were made during four ob- 
serving runs from 1995 to 2003 with the Swedish-ESO- 
Submillimetre Telescope SEST at the La Silla observatory, 
Chile. The SEST 3 and 2 mm (SESIS) dual SIS single side- 
band (SSB) receiver was used. The observations were made 
in the dual beam switching mode (beam throw 1L5 in az- 
imuth). The SEST half-power beam width (HPBW) and the 
main beam efficiency, rj MS , at frequencies 86 GHz, 115 GHz 
and 147 GHz are 57", 45", 34" and 0.75, 0.70, 0.66, respec- 
tively. The SEST high resolution 2000 channel acousto-optical 
spectrometer (bandwidth 86 MHz, channel width 43 kHz) was 



split into two halves to measure two receivers simultaneously. 
At the observed wavelengths, 2 mm and 1 mm, the 43 kHz 
channel width corresponds to approx. 0. 12 km s _1 and 0.08 km 
s _1 , respectively. The observed molecular transitions, their rest 
frequencies and the upper level energies are listed in TableQ] 

Calibration was achieved by the chopper wheel method. 
Pointing was checked regularly towards known circumstellar 
SiO masers. Pointing accuracy is estimated to be better than 
5". 

The 28 SiO (v = 0, J = 2 - 1) and (v = 0, J = 3 - 2) 
observations were carried out in October 1995 and April 1996 
and are described in detail in HLBZ98. 29 SiO(2 - 1) and (3 - 2) 
observations were carried out in October 1998. Typical values 
for the effective SSB system temperatures and the RMS noise 
of the spectra were 140 K and 170 K, and 0.016 K and 0.021 
K at frequencies 86 GHz and 130 GHz, respectively. Linear 
baselines were subtracted from the spectra. 

The CH3CCH (J = 5 K - A K ) and (J = 6 K - 5 K ) obser- 
vations were carried out in June 2003. At this time the elec- 
tronics of the SESIS 2 mm receiver had been decommissioned 
and only the 3 mm receiver was available. These were the last 
spectral line observations made with SEST. For each source 
typically six spectra, each with 2 minutes integration time, 
were obtained. The effective SSB system temperatures were 
130 K for the CH 3 CCH(5^ - 4 K ) transition and 250 K for the 
CH 3 CCH(6k - 5 K ). The RMS of the spectra were 0.017 K and 
0.024 K for the two transitions, respectively. Because of the 
better performance of the receiver at the lower frequency, the 
J = 5 k - 4-k transition was chosen for the survey. A third order 
baseline was subtracted from the spectra. Four ^-components 
(K = 0, 1, 2, 3) were detected towards all sources, and Gaussian 
profiles were fitted to the lines. The CLASS program of the 
GILDAS software package developed at IRAM was used for 
the reductions. In the fitting procedure the line separations be- 
tween the components were fixed and the line widths for all 
/^-components were assumed to be identical. 

2.2. Continuum 

The 1.2 mm continuum observations were carried out in June 
2003 with the 37 channel SEST imaging bolometer array, 
SIMBA. The SIMBA central frequency is 250 GHz and the 
bandwidth about 50 GHz. The HPBW of a single bolome- 
ter element is ~ 24" and the separation between elements 
on the sky is 44". The observations were conducted in stable 
weather. Frequent skydips were used to determine the atmo- 
spheric opacity and the values obtained varied between 0.27 
and 0.3. Pointing was controlled by observing sources with 
known accurate coordinates and is estimated to be better than 
5". Uranus was used for flux density calibration. 

The observations were done in the fastscanning mode with 
a scanning speed of 80" s _1 . The typical map consisted of 65 
scans of 700" in length in azimuth and spaced by 8" in ele- 
vation. Each source was observed at least twice. The two maps 
were done at different LST values to reduce the possibility of 
strong artefacts in the reduced maps. The data were reduced 
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using the MOPSI 1 software package according to guidelines 
in the SIMBA Observers Handbook (2002) and Chini et al. 
(2003). 

All the observed maps contained strong sources and there- 
fore, as suggested in the SEST handbook, no despiking was 
applied to the data. The data reduction deviated from the sug- 
gestions in the SEST handbook in the following details. The 
SEST dish shape (and the beam profile) degrades as a function 
of the elevation as the dish deforms because of its own weight. 
Gain-elevation correction is used to compensate the decreas- 
ing of the observed source peak intensity as a function of the 
telescope elevation. However, no gain-elevation correction to 
the present data was applied because aperture photometry was 
used to define the source flux. The source flux lost because of 
the decrease of the telescope gain in the source position is re- 
covered as the source is mapped. The telescope gain is unity 
at elevation of 74.6 degrees and decreases by 10 % at 45 de- 
grees. Most of the present observation were obtained at eleva- 
tions 50 degrees or higher where the said correction is small. 
Third order baseline was applied to the data before correlated 
noise removal after which all the bolometer elements with in- 
sufficient base range were masked. A first order baseline was 
applied to individual scan lines after this. The refinement of the 
SIMBA data reduction is discussed in detail in Haikala (2006, 
in preparation). 

3. Observational results 

The observed sources are listed in Table |2 The columns of 
this table are: (1) source name; (2) and (3) equatorial coor- 
dinates (J2000.0); (4) galactic coordinates; (5) distance; (6) 
galactocentric distance; (7) LSR velocity determined from the 
CH3CCH line; (8) notes on association with molecular masers 
or an UC Hn region. When no distance reference is given 
in Col. 5, the value is the kinematic distance calculated from 
the CH3CCH velocity using the rotation curve determined by 
IB rand & Blitz 19931 and Rq = 8.5 kpc (galactocentric distance 
of the Sun) and ©o = 220 km s _1 (circular velocity at a distance 
Ro). 

3.1. Dust emission 

The obtained SIMBA maps are presented in Fig. [2 The maps 
are plotted to the same angular and intensity scales. In Table 
[5] we list the coordinates (0:2000.0, 62000.0) an d the intensities 
(7™ ax ) at the dust emission maxima identified on the maps. The 
SEST HPBW at the frequency of the continuum observations 
is 24"which is less than half of the HPBW at the frequency of 
the observed SiO 7 = 2-1 and the CH3CCH lines. In order 
to estimate the average dust and H2 column densities within 
the beam used for the 3 mm line observations, the SIMBA 
maps were smoothed to the resolution of 57". The resulting 
smoothed surface brightness (7™°) towards the position of the 
line measurements is listed in Col. 2 of Table 0] The 1.2 mm 
flux density (S 1.2 mm) integrated over the source and the angu- 

1 MOPSI is a software package for infrared, millimeter and radio 
data reduction developed and constantly upgraded by R. Zylka. 



lar size (FWHM) of the core (0 S ) determined from the orig- 
inal map is listed in Cols. 3 and 4 of this Table. The angu- 
lar sizes were estimated from two-dimensional Gaussian fits to 
the surface brightness distribution (geometric mean of the ob- 
served major and minor axes of Gaussians). The source size 
was corrected for the broadening by 24"the telescope beam in 
accordance with the assumption of Gaussian distributions. The 
flux density was calculated in an 72"diameter circular aperture 
which correponds to three times the original beam size. 

3.2. SiO and CH 3 CCH lines 

The obtained SiO spectra are presented in Fig. 13 The lines are 
nearly Gaussian in the central part and have broad, asymmet- 
ric wings. The line parameters are listed in Table [6] In this 
table we give the minimum, maximum, and peak velocities 
(pmin, fmax an d u pe ak, respectively), and the integrated intensi- 
ties, J T^(v)dv, over the given velocity range. 

The spectra of the two SiO isotopologues were used to es- 
timate the optical thicknesses of the 7 = 2-1 and 7 = 3-2 
lines of 28 SiO, ti^\ and T3^2- For this purpose all spectra were 
resampled to the same LSR velocity grid with a channel width 
of 1 km s _1 . The smoothing of the spectra to this resolution was 
necessary to achieve the sufficient signal to noise ratio. Using 
these resampled spectra the optical thicknesses T2^i and T3^2 
could be derived in 1-13 velocity channels near the line peak 
for twelve sources. The terrestrial SiO isotopic abundance ra- 
tio is X tel -( 28 SiO/ 29 SiO) = 19.6. In interstellar medium values 
in the range 10-20 have been determined (Penzi as"l98U . Here 
we use the value X = 20 when determining the optical thick- 
nesses. The 28 SiO optical thicknesses T2^i and T3^2 are typi- 
cally between 1 and 3 in the line centres. 

The excitation temperatures, 7 ex , were estimated from the 
optical thickness ratio by using Eq. (2) of HLBZ98. These es- 
timates are based on the assumption that 7 ex (7 = 3-2) = 
7 ex (7 = 2-1). They do not depend on the beam filling as the 
optical thicknesses were determined using pairs of lines having 
similar frequencies. The derived values of 7 ex are close to 5 K 
for all sources and velocity channels. The weighted average of 
J ex and its standard deviation are 7 ex = 4.3 ±0.1 K (see Table 
®. 

From these results, the SiO column densities, A^(SiO), were 
estimated from the integrated intensities of the 29 SiO spectra 
by assuming optically thin emission and uniform excitation of 
the rotational lines with 7 ex = 5 K for all sources. 

The CH3CCH spectra are shown in Fig. [3] and the 
Gaussian parameters of the detected lines are given in Table 
[7] The line peak velocities and widths (FWHM) are listed 
in Cols. 2 and 3 of this table. The integrated intensities of 
the components K = 0,1,2 and 3 are given in Cols. 4 - 
7. The integrated intensities were used to derive the rota- 
tional temperatures, T Iot , and the CH3CCH column densities, 
Af(CH3CCH), with the population diagram method described 
and discussed in, e.g., lAskne et al. 1984| |Bergin et al. 1994| 
and |Goldsmith & Langer 1999] 
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Fig. 1. SIMBA images of the cores (1). The grey scales and contours represent the intensity of 1.2 mm dust emission. The 
contour levels are 0.25, 0.5, 1.0, 2.0, 4.0, 8.0 and 16.0 Jy/beam in all maps. The maps are plotted to the same angular scale. The 
target positions of the line observations are denoted with crosses. The filled circles indicate the locations of mid- and far-infrared 
sources. The scale bar on the left corresponds to 1 pc and the beam HPBW is shown in the lower right corner. The image borders 
have been masked out. 



4. Physical/chemical parameters of the cores 

4.1. Column densities and rotational temperatures 

In the case of uniform excitation and optically thin dipole tran- 
sition, the total column density of the emitting molecule is re- 
lated to the integrated line intensity by the formula 



3ke Q 1 e E -^Z(T ex ) 1 



2tt 2 /j 2 vS u i 



F(T be ) 
f(7-„) 



1 j T* A (v)dV: 



(1) 



where k is the Boltzmann constant, e () is the vacuum permit- 
tivity, fi is the permanent dipole moment, v and S u are the fre- 
quency of the transition and the line strength, E u is the upper 



state energy, Z is the rotational partition function, rj is the beam- 
source coupling efficiency, and the function F(T) is defined by 

1 



F(T) = 



(2) 



e hv/kT _ I ' 

For the rotational transition J u — > J u - 1 of a linear molecule, 
like SiO, S u = J„- We have assumed that r ex = 5 K for all 
sources (see Table |5ji and 77 = 77MB (0.75 at 86 GHz and 0.68 
at 129 GHz), which means that the source just fills the main 
beam. For the SiO isotopic abundance ratio we have used the 
value X = 20. 

The rotational temperatures, T mt , and the CH3CCH column 
densities, A^(CH3CCH), were derived by means of the popu- 
lation diagram method. This method is based on comparison 
of intensities of spectral lines which lie close in frequency al- 
though arising from rotational levels with different energies. In 
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the case of CH3CCH, the close lying transitions represent dif- 
ferent /^-components. 

Assuming uniform excitation and optically thin emission, 
the 'rotational diagram equation', i.e. the equation relating the 
integrated intensities of different Zf-components to the rota- 
tional temperature, T ro t, and the total column density, N tot , can 
be written as 



In 



= In 



27rV Mc 
3keo Zjd 



J_Eu 



(3) 



where g^ is the K degeneracy and statistical weight g; takes 
the spins of hydrogen nuclei into account. See, e.g., Appendix 
B of Anders on et al. 199 9 for derivation of Eq. 0. As above, 
we have again assumed that 77 = 77MB- 

Rotational diagrams for two sources are shown in Fig. |4] 
Straight lines were fitted to the data using a least-squares fit- 
ting technique. The data are consistent with a single T TOt . The 
goodness of the fit substantiates the assumptions that the gas is 
in LTE with T mt as 7\i„ and that the lines are optically thin. 

The H2 column density, NQI2), of each source was derived 
from dust continuum emission using the following equation: 



JV(H 2 ) = 



B v (Td)ymHKdRd 



(4) 



where /jj llst is the observed dust surface brightness (see Table 
0}, B Y {TS) is the Planck function for a blackbody of dust tem- 
perature 7d, // = 2.29 is the mean molecular weight account- 
ing for a 10 % contribution of helium (e.g.. IHill et al. 20 05 1. 
mh = 1.673534 • 10~ 27 kg is the mass of the hydrogen atom, k& 
is the mass absorption coefficient per unit mass of dust, and Rd 
is the dust-to-gas mass ratio. 

We assume that 7^ is equal to the gas kinetic temperature 
derived from CH3CCH observations. T^ n and 7^ are believed 
to be well coupled deep in the cloud, above a density of about 
10 4 cm~ 3 , due to collisions between gas and dust grains (see, 
e.g., Takahashietal. 19831 PeetaT 2004 1. Values of 0.1 m 2 



kg 1 at A — 1.2 mm (Ossen kopf & Henning 1994) and are 
adopted for and /?a, respectively. 

The rotational temperatures, the H2, SiO and CH3CCH 
column densities and the fractional SiO and CH3CCH abun- 
dances (x(x) = N(x)/N(H2)) are listed in Table [8] The av- 
erage value and the standard deviation of T rot is 32.7 + 3.6 
K. The N(U 2 ) values are found to be ~ 10 22 - 10 23 cnr 2 , 
and /V(SiO) and /V(CH 3 CCH) lie in the range 0.8 - 5.5 • 10 13 
cm" 2 and 2.0 - 15.0 • 10 14 cm~ 2 , respectively. The values of 
*(SiO) and ^(CH 3 CCH) are found to be ~ 3.0 + 1.7 ■ 10 10 and 
6.9 ±4.1 • 10~ 9 , respectively. 

4.2. Linear sizes, mass estimates and densities 

The linear sizes have been computed from the angular diame- 
ters listed in Table|4]using the distances listed in Table|2] 

The core masses have been estimated using the 1.2 mm 
continuum maps. We have also calculated the virial masses us- 
ing the velocity dispersions and kinetic temperatures from the 
CH3CCH data. The mass of a core, M cont , has been calculated 
assuming that the dust emission is optically thin and that the 
dust-to-gas ratio, Ra, and the dust absorption coefficient per 
unit mass, /Cd, are constant. These assumptions imply the fol- 
lowing equation: 



S v d 2 



BviT^Rt 



(5) 



where S v is the 1.2 mm continuum integrated flux within the 
selected aperture (see Tabled*. 

The virial masses, M v - U -, have been estimated by approx- 
imating the mass distribution by a homogenous, isothermal 
sphere without magnetic support and external pressure, using 
the formula 

M vir [M Q ] = (0.5(A i ; 1/2 [kms- 1 ]) 2 +0.01 T[K])d[kpc]® s ["] ,(6) 

where Ai>i /2 is the FWHM of the CH 3 CCH(5^ - 4 K ) line, T is 
the gas kinetic temperature derived from the CH3CCH obser- 
vations, d is the distance to the source, and S = 2R/d is the 
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Fig. 1. continued. 



angular diameter of the source as estimated from the SIMBA 
maps. For W28 A2(2) no virial mass estimate was obtained as 
its (deconvolved) angular size is less the half the beam size. 

The average H2 number densities, n(H2), were calculated 
using the masses, M cont , and radii, R, estimated from the dust 
continuum maps. The obtained radii, masses and average den- 
sities are listed in Table|5] The average radius of the sources is 
0.2 + 0.16 pc, the average value of the masses estimated from 
continuum emission is 1.8 + 2.1 • 10 3 M©, and the average den- 
sity is 4.5 + 6.5 • 10 6 cirT 3 . 

5. Discussion 

5.1. SiO and CH 3 CCH abundances 

The 1 .2 mm dust emission is likely to be dominated by the cool, 
T ~ 30 K, envelopes, although it has a contribution from hot 
cores (T > 100 K) around newly born massive stars. The dust 
temperatures derived by Fau ndez et al. 2004l towards several of 
our sources using SIMBA and IRAS 100 and 60 fim data are 
comparable to the CH3CCH rotational temperatures, r rot , de- 
rived here. Taken the uncertainties due to different beam sizes 



of the SIMBA and IRAS into account the agreement is reason- 
able and substantiates the assumption Ti ust « T mt used in the 
total NQI2) estimates. 



While it can be assumed that the CH3CCH lines and 1.2 
mm dust continuum emission trace to a large part the same 
material, this is not obvious for dust and SiO. All SiO lines 
observed in this survey show high-velocity wings, are single- 
peaked, and have their maxima near the systemic velocity of the 
cloud. These characteristics can be qualitatively explained with 
bow-shock models where the jet generating the shock is seen 
in a small angle, i.e. either head on or tail on. A large inclina- 
tion should result in a double peaked profile from a bow shock, 
and in a narrow line in the case where the emission originates 
in a turbulent wake behind a bow shock (HLBZ98, Fig. 12). In 
view of the distance to the sources and the angular resolution 
of the present observations, it is clear that the telescope beam 
encompasses entire star-forming regions with various gas com- 
ponents and young stars at different evolutionary stages. The 
line profiles may therefore have contribution from several out- 
flows with different inclinations. 
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In the light of previous SiO mapping observations it seems 
possible, however, that part of the low-velocity emission orig- 
inates in the quiescent envelope. This kind of component 
has been discovered in high-resolution SiO mappings of star- 
forming regions (Lefloch et al. 1998; Codella et al 1999; 
Shepherd et al. 2004; Fuente et al. 2005). In these stud- 
ies SiO emission characterized by narrow and broad lines 
have been found to come from separate regions. The SiO 
abundances derived for the low-velocity component are ~ 
10 -10 . This value is smaller than the abundances derived to- 
wards high-velocity shocks (k, 10~ 8 , e.g., Bachi ller et al. 199TI 
Gib b et al. 2004> . and larger than the upper limits derived in 
cold, dark clouds and PDRs (~ 10~ 12 , e.g.,|Ziurys et al. 1 989 
iMartrn-Pintado et al. 1 992*1 ISchilke et al. 20oH . 

In the present single-pointing observations typically half of 
the integrated SiO intensity comes from the velocity range with 
detectable CH3CCH, which is supposed to trace the quiescent 
gas component (see Table[6j. However, the same radial velocity 
does not necessarily mean spatial coincidence. Although SiO 



abundances derived in outflows are based on comparison with 
CO in the same velocity range, a similar spatial velocity cor- 
relation cannot be assumed near the cloud's systemic velocity. 
Depending on the structure and orientation of the shock, it is 
possible to find shock-produced SiO at low radial velocities. 
Therefore the column density of the low-velocity SiO gives 
only upper limit to the SiO column density in the quiesent en- 
velope, and the same is probably true for the fractional SiO 
abundance using the H2 column density determined by dust. 

A comparison of the SiO main beam brightness tempera- 
tures, Tmb, with the derived optical thicknesses, t, and excita- 
tion temperatures, T ex , suggests that the beam filling factors are 
close to unity, and the SiO emission regions are extended. We 
find that the upper limits for the fractional SiO abundance in 
the low-velocity gas are on the order of 10 -10 , i.e. similar to the 
SiO abundances derived previously for quiescent gas compo- 
nent in some star-forming regions. It should be noted, however, 
that the same result would be achieved if SiO is present only in 
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Fig. 1. continued. 



shocked layers filling 1% of the gas volume. This is not neces- 
sarily in contradiction with nearly uniform beam filling. 

The model for greatly enhanced SiO production 
in powerful shocks by silicate grain destruction and 
subsequent high-temperature gas-phase chemistry is 
well established (ISchTlke et al. 19971 ICaselli et al. 19971 
IPineau des Forets et al. 1997b . The possible occurrence of SiO 
in the quiescent gas component is not as well understood. In 
the postshock gas SiO should convert in a reaction with OH 
into Si0 2 (SiO + OH — > Si0 2 + H), which is eventually 
removed from the gas phase by accretion onto the dust grains. 
According to Cod ella et al. 1 999 the times needed to remove 
SiO from the gas phase and to slow down speeding SiO 
'blobs' are similar, ~ 10 4 yr, and this would support the 
outflow remnant scenario. In another, perhaps more feasible 
model high-velocity SiO emission comes from bow-shocks 
generated by high-velocity jets. In this case low-velocity SiO 
can originate in turbulent wakes behind bow-shocks (see 
Rag a and Cabrit 1993| and the model profiles in Fig. 12 of 
HLBZ98). 

It has been suggested that SiO can be produced also in 
warm, quiescent gas via neutral or ion-neutral reactions, possi- 
bly preceded by evaporation of Si-bearing molecules from the 
icy mantles of dust grains. The suggestion of Ziurys et al. 1989 
and |Langer & Glassgold 1990 that neutral reactions with acti- 
vation energies are important for the SiO production in warm 
gas, would imply a strong correlation between the SiO abun- 
dance and the average kinetic temperature. As shown in Fig. 
|5J bottom panel, even in the rather narrow temperature range 
covered, 25 - 39 K, a substantial change in the fractional 
SiO abundance would be expected if it were proportional to 
exp(-lllK/r) as suggested by Langer& Glassgold 1990 In 
the model of MacKay i |MacKay 1995||MacKay 1996) icysili- 



con is mainly in the form of SiH4, which can evaporate when 
the dust is warmed up. Also in this model one could expect the 
warmer cores have larger SiO abundances. 

The abundances derived for the low-velocity SiO do not 
show, however, any positive correlation with the average ki- 
netic temperature. Instead, they even seem to decrease slightly 
when the temperature rises (Fig.|5] bottom panel). 

Also the CH3CCH production can be enhanced by inten- 
sified desorption. CH3CCH is one of the lowest energy prod- 
ucts in reaction between the methylidine radical (CH) and 
ethylene (C2H4, Cano sa"et al. 1997l > which is formed on dust 
grains. The diagram presented in Fig. |5J top panel, suggests 
a positive correlation between the fractional CH3CCH abun- 
dance and 7\i n , which would conform with the desorption sce- 
nario. On the other hand, this tendency does not agree with 
the prediction of the pure gas-phase model of ILee et al. 1996 
(see also Ala koz et al. 20 02 ), which predict a strong anticorre- 
lation between the CH3CCH abundance and the kinetic tem- 
perature. The CH3CCH abundances derived here are simi- 
lar to those found in the Orion Ridge, M17 and Cepheus 
A ( |Ungerechts et al. 1997( |Bergin et al. 1997) , and substan- 
tially larger than those determined towards cold, starless cores 
JKontinen et al. 20001 IMarkwick et al. 20051 . 

These findings lend support to the possibility that indeed 
grain-surface reactions are important for the production of 
CH3CCH and its gas-phase abundance depends on temperature 
via the thermal desorption. On the other hand, no evidence for 
such behaviour is found in the case of SiO. This complies with 
the shock origin of SiO. The diminishing tendency towards 
warmer cores suggested by Fig.|5]may reflect an evolutionary 
effect. 

According to the picture presented by IFuente et al. 2005 
based on a study of the environments of two young stars in 
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Fig. 2. The v - 0, 7 = 3- 2 and y = 0, 7 = 2-1 lines of 28 SiO and 29 SiO towards each observed source. The weaker 29 SiO 
spectra are multiplied by 2. 
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Fig. 3. Spectra of the CH3CCH J = 5k - 4-k and J - 6k - 5k rotational transitions towards each observed source. The X-axis 
represents the LSR velocity of the 5 1 - 4j line (62 - 52 in the case of 6k - 5k) and the intensity scale is antenna temperature, 7^, 
in Kelvins. The fairly low line intensities, T* A < 1 K, suggest a low optical thickness. 



the region of the nebula NGC 7129, protostellar envelopes are 
dispersed and warmed up during the early stellar evolution, at 
the same time as shocks associated with outflows get less en- 
ergetic. For SiO this means decreasing abundance with time. 
The abundance variation from core to core is not very large in 
our sample. This is probably related to the fact that all the stud- 
ied objects are associated with powerful masers and/or UC H 11 
regions, and therefore the dynamical and chemical age vari- 
ations are likely to be rather limited ( |Garay & Lizano 1 999 
|Bergin et al. 1997) . 

5.2. Dynamical states 

Most of the cores studied here appear to be gravitationally 
bound which is indicated by the fact that M cont ^ M vlr . 
The M cont /M v i r ratios range from 0.7 to 6.8, the average be- 
ing 2.5. For five cores M cont is substantially larger than M v h-, 
whereas the rest appear to be close to the virial equilibrium 
(M v h- ~ M cont , within a factor of two). The virial masses, M v - U -, 
were calculated assuming homogenous density distributions. 
A more realistic distribution described by a powerlaw of the 
type n(r) oc r~ p , where p is typically ~ 1.5 - 2.5, would make 
the virial mass estimates smaller (see, e.g,. IFontani et al. 2005> . 
On the other hand, the effect of the density gradient can be 
compensated by the increase of the kinetic temperature to- 
wards the core centre. All the cores studied here contain newly 
born massive stars which are known to give rise to expand- 
ing UC Hn regions and hot cores. Magnetic fields are likely 
to provide additional support against the gravitational collapse 
(e.g., |Garay & Lizano 19 99 Font ani et al. 20021 1. For example, 



Sarma et al. (2000) detected magnetic fields toward sources A, 
E, and D in the NGC 6334 complex, which are of the same 
order as the critical fields required to support the cores against 
gravitational collapse. 

5.3. Comments on individual sources 

Ten of our fourteen targets have been included in re- 
cently published, extensive millimetre continuum surveys with 
SIMBA. Fau ndez et al. 20041 have mapped about 150 mas- 
sive cores around IRAS point sources with FIR colors typ- 
ical of UC Hn regions, and discuss the general characteris- 
tics of these objects. They present the maps of IRAS 12326- 
6245, G326.64+0.61 (IRAS 15048-5356), OH328. 81+0.63 
(IRAS 15520-5234), G330.95-0. 19 (IRAS 16060-5146), IRAS 
16562-3959, G345.00-0.23 (IRAS 17016-4124), NGC 6334F 
(IRAS 17175- 3544), and G35 1.77-0.54 (IRAS 17233-3606). 
The survey of IHill et al. 2005 covering about 130 regions in- 
clude IRAS 12326-6245 (G301. 14-0.2), G330.95-0.18, W28 
A2 (G5.89-0.39), and W31 (2) (G10.62-0.38). 

The structure of the sources in Fig. ^ agrees well with 
those presented in Faun dez et al. 2 004 and Hill et al. 20051 The 
peak fluxes in IHill et al. 2005l for the sources in common with 
this work agree also well. The integrated flux densities for the 
sources are however higher in lHill et al. 200 5 and especially in 
Fau ndez et al. 2004l than in the present paper. Hill et al. calcu- 
lated the integrated flux density inside a polygon containing all 
the emission around the source whereas in the present paper 
the integration was done inside a circle corresponding to three 
times the telescope beam width. Contrary to the present pa- 
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Fig. 4. Rotational diagrams of the CH3CCH spectra for two of the observed sources. X-axis plots the upper transition state 
energies divided by the Boltzmann constant, -f-, and Y-axis plots the lefthand side of Eq. 0. The lines represent linear least- 
squares fits to the data. The error bars correspond to the errors of integrated intensity according to the law of propagation of 
errors. The rotational temperature, r rot , and the total CH3CCH column density, N, are given in the top right of each panel. 



per IHill et al. 2 005 also applied the gain- elevation correction 
to the SIMBA data. The different integration aperture and the 
gain-elevation correction can well explain the Hill et al. (2005) 
~ 30 - 50 % higher flux densities. Faunde z et al. 200 4 give no 
details on how the integrated fluxes were calculated so it can 
not be decided if the higher values are due to e.g. the integra- 
tion aperture. 

The kinetic temperatures from CH3CCH are lower than 
the dust temperatures derived by Fau ndez et al. 20041 from the 
spectral energy distributions (SEDs) using the flux densities 
at 1.2 mm and in the four IRAS bands. However, as our 1.2 
mm flux densities are, in general, lower than those derived 
by Faundez et al., we end up with comparable masses. The 
masses derived by IHill et al. 20051 are on the average larger 
than those derived by us, mainly because they have obtained 
slightly higher flux densities and assumed that = 20 K. 
Finally, the mass estimates are proportional to the distance 
squared, and there are some differences in the derived kine- 
matic distances depending on the Galactic rotation curve and 
y LSR used. 

Next we discuss six of our sources, IRAS 12326-6245, 
NGC 6334, G35 1.77-0.54, W28 A2(2), W31 and W33 in 
more detail. In general, the sources are typical high-mass star- 



forming cores associated with several molecular masers, UC 
Hn regions and infrared sources (see Tableland Fig.Q. 

5.3.1. IRAS 12326-6245 

IRAS 12326-6245 is a luminous FIR source (L bol = 3.8 • 10 5 
Lq) which is located at a kinematically estimated distance of 
4.4 kpc JZinchenko et al. 19951 IDsterloh et al. 1997b . This ob- 
ject is associated with two UC Hn regions, both of which 
are associated with mid-infrared sources Penning et al. 2000) . 
Also several molecular masers including H2O, OH and CH3OH 
have been identified at the position of IRAS 12326-6245 (see 
Pfenning et al. 2000| for references). 

The molecular line maps (see Henning et al. 2000 and ref- 
erences therein) shows that one of the most energetic and mas- 
sive bipolar molecular outflow (mass outflow rate M ~ 0.02 
Mq yr _I ) in the southern sky is originated close to IRAS 
12326-6245. 

The core 1.3 mm flux density and the gas mass obtained by 
Henning et al. (12.0 Jy and 2400 Mq, respectively) are very 
close to the values in the present study (12.8 Jy and 2200 Mq, 
respectively). 
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Fig. 5. Fractional CH3CCH and SiO abundances as functions 
of the rotational temperature, T rot , derived from CH3CCH. The 
solid lines represent linear least-squares fits to the data. The 
dashed line in the bottom panel represents a change in the frac- 
tional SiO abundance proportional to exp(-lllK/r) as sug- 
gested by |Langer & Glassgold 1990] 

5.3.2. NGC 6334 

The star-forming molecular ridge associated with the giant 
Hn region NGC 6334 is the best studied of our targets. For 
recent, comprehensive molecular line and continuum studies 
see Kraemer et al. 1999, Kraemer & Jackson 1999, Jackson 
& Kraemer 1999, Sandell 2000, McCutcheon et al. 2000 and 
Brooks & Whiteoak 2001. These authors summarize the ac- 
cumulated knowledge of the region, and include keys of the 
nomenclature. Sandell (2000) mapped the northern part of the 
molecular ridge with UKT 1 4 bolometer at the JCMT using four 
millimetre and submm wavelengths in the range 1.1-0.35 mm. 

The six brightest 1.2 mm peaks listed in Table [3] can 
be associated with star-forming cores detected as far-infrared 
sources or UC Hn regions as follows: 1) FIR-I/UC Hn 
F; 2) I(N); 3) FIR-IV/UC Hn A; 4) FIR-V; 5) FIR- 
III/UC Hn C; 6) FIR-II/UC Hn D JMcBreen et al. 19791 



Rodriguez et al. 1982| IGezari 1982V The 1.2 mm peak No. 
7 is located on the southeastern side of radio shell be- 
tween FIR-IV and V, in the region of the PDR G35 1.2+0.70 



(Jackson & Kraem e~1999> . It is associated with a CO clump 
in the survey of Krae mer & Jackson 19991 (their Fig. 5). The 
dust emission peak No. 8 on the northwestern side of the 
dense molecular ridge is associated with the object called 'dust- 
cloudlet' bv lSandell 20001 

Assuming a uniform dust temperature of 30 K, we obtain 
from the SIMBA map a mass of 1 1 400 Mq for the entire ridge. 
The share of its northern end with the cores I, I(N) and II is 
5400 Mq, and for the regions III, IV, and V we get 900 Mq, 
2200 Mq, and 1300 Mq, respectively. 

The main features of the 1.1 mm map of Sandell 2000 can 
be recognized on the present SIMBA map, despite a lower S/N 
available here. A careful inspection of the image reveals fila- 
ments reaching out from the dense ridge, bearing much like- 
ness to the structures associated with OMC-1 discovered by 
|Johnstone & Bally 1999|onaJCM T/SCUBA map. The sugges- 
tion of I Johnstone" & Bally 1999| that these features could rep- 
resent the cavity walls of past outflows is appropriate also in 
this case. Filamentary structures are rather common features 
in GMCs forming high mass clusters, so there is no reason to 
consider them as artefacs of the data reduction. 

The 'linear filament' discussed by ISandeTHZ OOO continues 
to the north and southwest and in fact the ridge is arc-shaped. 
In Fig. [6] an arc of a circle of radius 27' (13 pc) and centred 
at R.A. nMS^OiO, Dec. -35°35'30" is superposed on the 
SIMBA map. The alignment evident in the figure suggests that 
the ridge has been formed by an expanding Hn region with 
its centre near the given position. This H 11 region has proba- 
bly been ionized by one or more of the matured OB stars in 
the region. The closest O-type star on the western side of the 
ridge is HD 319 699 (07, distance 1.6 kpc). The UV radiation 
from this star is likely to dominate the excitation of the ex- 
tended emission nebula in this direction. It seems, however, 
that the star has not been in the centre of the expansion, be- 
cause it lies about 11' southeast of the centre of the circle, 
and its proper motion vector (although very uncertain) points 
northwest (The Tycho-2 Catalogue; |H0g et al. 20001 . The O- 
stars lying southeast and south of the ridge, HD 319702 (09, 
d = 1.4 kpc), HD 319 703 (06, d = 2.0 kpc), and HD 156738 
(07, d = 1.3 kpc) are surrounded by separate, roundish ioniza- 
tion regions which probably have not contributed to shaping the 
ridge (see, e.g., Fig. 5 of Kraemer & Jackson 1999 or SkyView, 
http://skyview.gsfc.nasa.g0v/1. 

5.3.3. G351 .77-0.54 

G35 1.77-0. 54 is the brightest SiO line source of the sample 
(see Fig. |2j. It is associated with the luminous steep-spectrum 
far-infared source IRAS 17233-3606, the strongest known 
ground-state OH maser source at 1665 MHz, other OH masers, 
and H 2 Q and CH 3 OH masers (e.g., ICaswell & Vaile 19951 
IMacLeod et al. 19981 IVal'tts et al. 20001 . The kinematic dis- 
tance derived using the LSR velocity from CH3CCH is 700 
pc, which would make this source the nearest in our sam- 
ple. The distance is smaller than the usually adopted value 
(~ 1.5 - 2.2 kpc), but it is close to the distance derived by 
IMacLeod et al. 19981 (1.0 kpc) from the LSR velocity of the 
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Fig. 6. The SIMBA dust continuum map of the star forming 
ridge in the region of NGC 6334. The ridge seems to follow an 
arc of a circle described in the text. The locations of two nearby 
O-type stars are indicated. 

H2CO 1 10- 1 1 1 absorption. The neighbourhood of the source di- 
rection to that of the Galactic centre causes a large uncertainty 
to the kinematic distance estimate. Nevertheless, the proximity 
is supported by the intensity of the SiO emission, suggesting a 
large beam-filling factor. 

5.3.4. W28 A2 

G5. 89-0.39 (also known as W28 A2) is one of the best-studied 
UC H 11 regions, but probably it does not belong to the best un- 
derstood. W28 A2 is associated with a highly energetic outflow, 
one of the most powerful ones in the Galaxy. Actually, in W28 
A2 there might be multiple outflows and driving sources. 

A bipolar outflow in the east-west direction with a mass 
of ~ 70 M© was identified by Harvey & Forveille (1988) 
in the J = 1-0 lines of CO, 13 CO, and C ls O toward this 
source. C 34 S observations shows a north-south oriented outflow 
dCesaroni et al. 199TV There is growing evidence that the 05 
star in G5.89 produced the N-S molecular outflow and hence 
is forming through accretion process Csee lArce et al. 2006l and 
references therein). The SiO (v — 0, J — 1-0) observations 
made with the VLA by Acord et al. (1997) shows a northeast- 
southwest bipolarity, which is confirmed by 1.3 mm contin- 
uum and several molecular line observations by Sollins et al. 
(2004). The 1.3 mm continuum source peaks in the centre of 
the Hn region ( S ollins et al. 20041 . Moreover, NIR observa- 
tions of the H2 v = 1-05(1) rovibrational line and the 
Bry emission show evidence of three outflows with distinct 
orientations and driving sources (Puga et al. 2006 1. Puga et al. 
suggest that the northwest-southeast oriented outflow is pos- 
sibly connected to the H2 knots seen in that direction. One 



of them ('knot B'), NW of W28 A2, can be identified in our 
SIMBA map. The two other features seen in our map, NW 
of knot B, at a 2 ooo = 18 h 00 m 15 s , 5 2 ooo = -24°01'30", and 
c*2000 = 18 h 00 m 08 s , 5 2 ooo = -24°00'30", could be related to 
NW-SE oriented outflow. 

Thus it is unlikely that there exist a single central high-mass 
star in W28 A2 as was initially thought to be the case due to 
the symmetric morphology of H n region. 

The kinematic distance derived from Hi observations to- 
wards the supernova remnant W28 by Velazquez et al. (2002) 
is 1 .9 kpc. Adopting this value instead of 2.7 kpc derived by us, 
the mass of the cloud would decrease from 970Mq to 480 Mq. 

Thompson & Macdonald (1999) found rotational tempera- 
ture of ~ 70 K from CH3CCH J = 20 - 19 and J = 21 - 20 ob- 
servations and Gomez et al. (1991) derived the temperature of 
~ 90 K from NH3(2, 2) and (3, 3) inversion transitions for the 
molecular envelope surrounding the UC H 11 region. We have 
observed much lower CH3CCH transition J = 5-4, which trace 
the gas in outer layers and thus giving lower rotational tempera- 
ture (34.6 K). The CH3CCH abundance estimate of 2.4- 10 7 by 
Thompson & Macdonald is an order of magnitude higher than 
our value of 1.5 • 10~ 8 . NH 3 (3, 3), (4, 4), and (5, 5) observations 
by Acord et al. (1997) imply the kinetic temperature of around 
100 K for the outflowing gas and they found the SiO abundance 
of 1 - 3 • 10~ 9 in this gas, which is also an order of magnitude 
higher than the abundance estimate in the present study in the 
position of W28 A2 (2). These results support the idea that the 
CH3CCH and SiO abundances are are clearly larger in hot core 
than in the warm envelope. 

5.3.5. W31 

W31 is one of the largest H 11 complexes in our Galaxy. G10.6- 
0.4 lying in the centre of our SIMBA map is one of the main 
extended Hn regions in this complex along with G 10.2-0. 3 
and G10. 3-0.1 (see, e.g., |Corbel and Eikenberry 2004| and ref- 
erences therein). NH3(1, 1) and (3, 3) observations show that 
the molecular core in G 10.6-0.4 is rotating and collapsing in- 
ward toward the UC H 11 region where molecular masers are 
seen distributed in the following way: H2O and CH3OH masers 
linearly in the plane of rotation, and OH masers along the rota- 
tion axis (see, e.g., Sollins & Ho 2005 and references therein). 
Further evidence for collapse comes from two-peaked profile of 
CO(4-3) line observed with the APEX ( |Wyrowski et al. 2006) . 
Observations made with the VLA by Sollins & Ho (2005) sug- 
gest that there is a rotating molecular accretion flow which flat- 
tens but does not form an accretion disk. Velocities of molec- 
ular lines observed from the surrounding neutral gas in G10.6- 
0.4 suggest that accretion flow passes through the H 11 region 
boundary and continues inward as an ionized flow ( Keto 2002 1. 
The associated infrared source (see Fig.0, IRAS 18075-1956, 
has a luminosity of 9.2 • 10 5 Lq (ICasoli et al. 1986i . 

The distance of W3 1 is not well known. However, spec- 
tral types of O stars identified in W3 1 and the NIR photome- 
try made by Blum et al. (2001) show that the distance to W31 
must be ^ 3.7 kpc. Corbel and Eikenberry (2004) suggest that 
G10.6-0.4 is located on the -30 km s _1 spiral arm at a distance 
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from the Sun of 4.5+0.6 kpc, which is consistent with the spec- 
trophotometric distance obtained by Blum et al. We have used 
the value of 3. 1 kpc which was derived by Blum et al. assuming 
that O star spectra are consistent with ZAMS. Note that Sollins 
& Ho (2005) and Sollins et al. (2005) adopt the kinematic dis- 
tance of 6.0 kpc from lDownes et al. 19801 

5.3.6. W33 

The W33 complex is heavily obscured in the visual. It 
contains a compact radio core G12.8-0.2. The FIR sources 
W33 A (G12.91-0.26) and W33 B (G12.70-0.17) associated 
with OH/H2O masers lie on the opposite sides of G12.8-0.2 
HCapps et al. 1978llStier et al. 19841 . 

IR and radio continuum observations have revealed a 
cluster of compact sources in the region of G12.8-0.2 
Pyck & Simon 1977| IHaschick & Ho 1983i . The total IR lu- 
minosity of G12.81-0.19 assuming a distance of 3.7 kpc is 
1 .4 ■ 10 6 L© JStieretal."l 984 1. According to the results of Stier 
et al. the minimum mass of the star associated with the radio 
core is 20 Mq. G12.8-0.2 is considered as a good candidate for 
a very young OB cluster (see, e.g. JBeck et al. 19 98 ). 

6. Summary and conclusions 

Spectral line observations at 3 and 2 mm, and dust contin- 
uum observations at 1.2 mm performed with the SEST were 
used to derive physical characteristics and fractional SiO and 
CH3CCH abundances in 15 high-mass star-forming cores. The 
sample represents typical GMC cores with the following aver- 
age properties < 7\ in >= 33 K, < n(H 2 ) >= 4.5 • 10 6 crrT 3 , and 
< M >= 1800 M©. While most cores seem to be near the virial 
equilibrium, for five of them the gravitational potential energy 
is large compared with the kinetic energy. 

The fractional SiO abundances determined from the ve- 
locity range with detectable CH3CCH emission were found 
to be ~ 1 - 7 • 10~ 10 . CH3CCH abundances were found to 
be 3 - 18 ■ 10~ 9 . In consistence with earlier results towards 
GMC cores (e.g., |Ungerechts et al. 1997| |Berginetal. 1997| > 
the CH3CCH abundance is found to be clearly larger than in the 
dense cores of dark clouds. It also shows a slight increase as a 
function of 7\; n . As this tendency does not agree with the pre- 
dictions of gas-phase chemistry models (e.g. lLee etal. 19 96 1, 
a possible explanation is that warm conditions lead to an inten- 
sified desorption of the precursor molecules of CH3CCH from 
the icy mantles of dust grains. This suggestion is supported by 
the chemistry model of Canosa ~tal. 19971 

The SiO abundances are midway between the upper limit 
from dark clouds and PDRs (~ 10 ) and the values de- 
rived for powerful shocks (~ 10~ 8 ). The abundance seem 
to decrease with rising temperature. This finding contra- 
dicts with the models where SiO production is dominated by 
neutral reactions with activation energies jZiurys et al. 1 989 
|Langer & Glassgold 1990) , and with models where desorption 
of Si-containing species from grain mantles are significant for 
the SiO production ( Ma c Kay 1996). As sug gested by the recent 
observational results of Fue nte et al. 20051 warmer cores repre- 
sent more evolved objects where highly energetic protostellar 



outlows releasing SiO into the gas phase are less frequent; this 
together with the rapid post-shock processing have lead to a 
diminished SiO abundance. 
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Table 1. The rest frequencies and the upper level energies of 
the observed lines. 



Table 3. Dust emission maxima indentified on the original 
maps (24" beam). 



Molecule 


Transition 


v [MHz] 


E a /k [K] 


2y SiO 


J = 


2 - 


- 1 


85759.0000 


6.195 


/s SiO 


J = 


2 - 


- 1 
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29 SiO 
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3 ■ 


-2 


128636.7064 


12.390 


28 SiO 
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3 ■ 


-2 
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12.504 


CH 3 CCH 


Jk = 


5 


-4 


85457.2720 
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CH 3 CCH 


Jk = 


5, 


-4, 


85455.6220 


19.505 


CH 3 CCH 


Jk = 


5 2 


_4 2 
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41.107 


CH 3 CCH 


Jk = 


5 3 


_4 3 
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CH 3 CCH 
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6 
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CH 3 CCH 
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6, 


-5, 
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CH 3 CCH 


Jk = 


62 


-5 2 


102540.1447 


46.029 


CH 3 CCH 


Jk = 


63 


-5 3 


102530.3487 


82.033 



Source 


peak No. 


«2000.0 


(52000.0 


I v [MJy/sterad] 
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1 
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2 
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Table 2. Source list. J2000.0 equatorial coordinates (a, 6), galactic coordinates (/, b), distance (d), galactocentric distance (/?gcX 
LSR velocity of the CH3CCH line (olsr) an d notes on associated masers and UC H 11 regions. 



Source 


ff2000.0 


#2000.0 


l,b 


d 


Rgc 


fLSR 


notes 




[h:m:s] 


[<•:':"] 


[°] 


[kpc] 


[kpc] 


[km s" 1 ] 
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H 2 F 


OH328.81+0.63 
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16:00:30.7 


-53:12:34 


329.03-0.20 


3.0 
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OH7H 2 O j /CH 3 OH g 


G330.95-0.19 


16:09:52.2 
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6.1 


4.3 


-100.0 


on H /n 2 o F 


G345.01 + 1.8N 


16:56:47.2 


-40:14:09 
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1.7 
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-13.8 
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IRAS 16562-3959 
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-40:03:42 
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-12.2 
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-27.7 


OH*/H 2 / 7CH 3 OH G 


NGC 6334 FIR-V 
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oh*7h 2 o f 


NGC 6334F 
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17:26:42.6 


-36:09:17 


351.77-0.54 


0.7 


7.8 


-3.2 


OH* r /H 2 F /CH 3 OH G /UC H 11" 


G353.41-0.36 


17:30:26.0 
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353.41-0.36 


3.5 
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-16.8 


OH M -*7H 2 0'7CH 3 OH , 7UC Hn" 


W28 A2(2) 


18:00:30.4 


-24:04:00 
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2.7 


2.4 


6.1 


OH / 7H 2 O jv /UCHii" 


W31 (2) 


18:10:28.7 


-19:55:50 


10.62-0.38 


3.1'' 


5.5 


-2.8 


OH / 7H 2 O a, /UCHii" 


W33 CONT 


18:14:13.6 


-17:55:25 


12.81-0.20 


3.7 


4.9 


34.8 


ftO" 



References: A ICohen et al. 19881 a ICaswell et al. 19891 c ICaswell et al. 19951 "Bronfman et al. 19961 £ IZinchenko et al. 19951 F The c atalog 

of non-stellar H 2 Q/OH masers JBraz & Epchtein 19831; G ICaswell 20001 g ICaswell 19981 7 ICaswell et al. 19951 ■' IScalise et al. 19891 K 
|Caswell & Haynes 1983| L |Moran & Rodriguez 1980|IGaume & Mutel 19871 M ICaswell et al. 19811 N Arcetri Atlas of galactic H 2 Q m asers 

<Comoretto et al. 1 990 Palag i et al. 1993|IBrand et al. 19941 . Distance references: " IZinchenko et al. 1995llOsterloh et al. 19971 * 
INeckel 19781 IBlum et al. 20011 When noreference is given, the distance is a kinematic distance calculated from the CH 3 CCH velocities. 



Table 4. Intensities in selected positions on the maps smoothed 
to 57", the total 1.2-mm flux density and the angular size 
(FWHM) of the source determined from the original maps 
(24" beam). 



Source I v [MJy/sterad] S h2mm [Jy] S ["] 



Table 5. Optical thicknesses of the 28 SiO(2- 1) and (3-2) lines 
and their excitation temperatures. In case the optical thickness 
ratios have been determined in several velocity channels (col- 
umn 5), the tabulated values are the averages over these chan- 
nels. 



Source 


T2- 


.1 


T5- 


•2 


T ex [K] 


chan 


G326.64+0.61 


3.0 


+ 


0.4 
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+ 


0.7 
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+ 


1.7 


2 


OH 328.81+0.63 


2.1 


+ 


0.2 


2.1 


+ 


0.2 


4.8 


+ 


0.4 


10 


IRAS 15566-5304 


3.8 


± 


0.4 


2.6 


+ 


0.6 


4.7 


+ 


0.6 


5 


G345.01 + 1.8N 


1.5 


± 


0.3 


1.3 


+ 


0.4 


5.3 


+ 


2.7 


1 


G345.00-0.23 


2.2 


± 


0.3 


1.5 


+ 


0.3 


4.6 


+ 


0.7 


6 


NGC 6334 FIR-V 


1.7 


± 


0.4 


1.1 


+ 


0.3 


4.1 


+ 


1.6 


1 


NGC 6334F 


2.3 


+ 


0.3 


1.2 


+ 


0.4 


3.4 


+ 


0.7 


3 


G35 1.78-0.54 


1.9 


± 


0.2 


1.7 


+ 


0.1 


4.2 


+ 


0.2 


13 


G353.41-0.36 


2.2 


± 


0.2 


2.0 


+ 


0.3 


6.1 


+ 


1.1 


5 


W28 A2(2) 


2.4 


± 


0.6 


2.0 


+ 


0.5 


5.2 


+ 


1.7 


3 


W31 (2) 


2.9 


± 


0.4 


2.9 


+ 


0.4 


4.0 


+ 


0.4 


6 


W33 CONT 


2.7 


± 


0.7 


3.0 


+ 


0.6 


7.2 


+ 


5.1 


1 



IRAS 12326-6245 


134 


12.8 


13 


G326.64+0.61 


98 


10.9 


27 


OH328.81+0.63 


136 


14.6 


22 


IRAS 15566-5304 


74 


7.8 


34 


G330.95-0.19 


154 


22.8 


13 


G345.01 + 1.8N 


154 


16.4 


47 


IRAS 16562-3959 


185 


19.1 


36 


G345. 00-0.23 


110 


11.4 


19 


NGC 6334 FIR-V 


189 


19.4 


39 


NGC 6334F 


360 


34.7 


16 


G351.77-0.54 


289 


27.5 


17 


G353.41-0.36 


197 


22.6 


33 


W28 A2(2) 


163 


15.5 


15 


W31 (2) 


198 


19.1 


89 


W33 CONT 


395 


41.1 


32 
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Table 6. 29 SiO integrated intensities for the line centres, i.e., for the velocity intervals with detectable CH 3 CCH emission, and 
for the whole line with wings. 









Line centre 








Whole line 






Source 


Line 




fmin 


^max 


fpeak 


J T* A (v)dv | cenlrc 


^min 


^max 


f T* A (v)dv | who l e 








[kms" 1 ] 


[km s" 1 ] 


[km s" 1 ] 


[Kkms" 1 ] 


[km s" 1 ] 


[krns" 1 ] 


[Kkm s~'] 


IRAS 12326-6245 


7 = 2- 


1 


-41.8 


-37.2 


-38.4 


0.10(0.01) 


-49.3 


-35.0 


0.20(0.02) 




7 = 3- 


2 






-40.0 


0.23(0.02) 


-48.4 


-32.0 


0.40(0.03) 


G326.64+0.61 


7 = 2- 


1 


-42.2 


-37.7 


-40.7 


0.26(0.01) 


-46.1 


-33.6 


0.47(0.02) 




7 = 3- 


2 






-41.0 


0.17(0.02) 


-45.5 


-33.6 


0.31(0.03) 


OH328.8 1+0.63 


7 = 2- 


1 


-44.6 


-38.3 


-41.9 


0.51(0.01) 


-55.3 


-24.9 


0.99(0.03) 




7 = 3- 


2 






-42.5 


0.60(0.02) 


-49.6 


-31.4 


1.07(0.03) 


IRAS 15566-5304 


7 = 2- 


1 


-46.2 


-41.0 


-45.8 


0.41(0.01) 


-53.1 


-33.4 


0.85(0.03) 




7 = 3- 


2 






-44.6 


0.26(0.02) 


-50.1 


-36.7 


0.50(0.03) 


G330.95-0.19 


7 = 2- 


1 


-92.7 


-88.7 


-91.6 


0.16(0.01) 


-99.6 


-80.5 


0.37(0.03) 




7 = 3- 


2 






-88.4 


0.13(0.02) 


-97.5 


-85.0 


0.20(0.03) 


G345.01 + 1.8N 


7 = 2- 


1 


-16.8 


-10.6 


-16.4 


0.31(0.02) 


-23.5 


-6.5 


0.50(0.03) 




7 = 3- 


2 






-17.3 


0.24(0.02) 


-21.1 


-7.4 


0.39(0.03) 


IRAS 16562-3959 


7 = 2- 


1 


-14.9 


-8.7 


-12.2 


0.21(0.02) 


-19.1 


-0.7 


0.35(0.03) 




7 = 3- 


2 






-11.3 


0.27(0.02) 


-21.8 


-2.4 


0.62(0.03) 


G345.00-0.23 


7 = 2- 


1 


-29.7 


-22.4 


-28.2 


0.43(0.02) 


-40.4 


-13.9 


0.88(0.03) 




7 = 3- 


2 






-26.9 


0.48(0.02) 


-40.7 


-11.2 


1.09(0.04) 


NGC 6334 FIR-V 


7 = 2- 


1 


-9.3 


-2.9 


-6.7 


0.23(0.02) 


-20.7 


-3.7 


0.51(0.03) 




7 = 3- 


2 






-7.0 


0.29(0.02) 


-17.7 


1.4 


0.56(0.03) 


NGC 6334F 


7 = 2- 


1 


-8.5 


-2.7 


-8.8 


0.22(0.01) 


-11.1 


-2.5 


0.34(0.01) 




7 = 3- 


2 






-10.1 


0.20(0.01) 


-13.8 


-0.1 


0.37(0.02) 


G351.77-0.54 


7 = 2- 


1 


-5.1 


0.8 


-3.9 


0.72(0.01) 


-19.0 


11.0 


1.85(0.03) 




7 = 3- 


2 






-3.7 


0.91(0.02) 


-17.8 


11.0 


2.29(0.03) 


G353.41-0.36 


7 = 2- 


1 


-19.7 


-13.6 


-16.7 


0.39(0.01) 


-24.7 


-6.2 


0.66(0.02) 




7 = 3- 


2 






-17.1 


0.30(0.01) 


-21.7 


-9.2 


0.42(0.02) 


W28 A2(2) 


7 = 2- 


1 


6.5 


12.8 


10.2 


0.23(0.02) 


-2.5 


26.6 


0.69(0.04) 




7 = 3- 


2 






11.9 


0.35(0.02) 


0.1 


37.9 


1.34(0.05) 


W31 (2) 


7 = 2- 


1 


-5.2 


1.6 


-4.5 


0.39(0.02) 


-12.8 


7.2 


0.68(0.03) 




7 = 3- 


2 






-2.1 


0.41(0.02) 


-12.5 


5.7 


0.67(0.03) 


W33 CONT 


7 = 2- 


1 


31.9 


39.2 


35.0 


0.24(0.02) 


28.9 


48.0 


0.33(0.03) 




7 = 3- 


2 






35.4 


0.28(0.02) 


27.1 


47.7 


0.48(0.03) 



Table 7. CH3CCH line parameters obtained from Gaussian fits. 



Source u LSR [kms 4 ] Aui/2 [km s" 1 ] J T A (v)dv [K-kms" 1 ] 









K=0 


K=l 


K=2 


K=3 


IRAS 12326-6245 


-39.54(0.04) 


3.68(0.02) 


0.69(0.02) 


0.58(0.02) 


0.29(0.02) 


0.13(0.02) 


G326.64+0.61 


-39.76(0.02) 


3.35(0.02) 


1.41(0.02) 


1.10(0.02) 


0.47(0.02) 


0.19(0.02) 


OH328. 81+0.63 


-41.72(0.01) 


3.07(0.02) 


1.89(0.02) 


1.57(0.02) 


0.63(0.01) 


0.30(0.01) 


IRAS 15566-5304 


-43.64(0.03) 


3.73(0.03) 


0.95(0.03) 


0.77(0.03) 


0.31(0.03) 


0.11(0.02) 


G330.95-0.19 


-90.95(0.05) 


5.42(0.06) 


1.31(0.02) 


1.00(0.02) 


0.54(0.02) 


0.21(0.02) 


G345.01 + 1.8N 


-13.78(0.02) 


3.60(0.02) 


1.23(0.02) 


0.96(0.02) 


0.48(0.01) 


0.24(0.01) 


IRAS 16562-3959 


-12.17(0.01) 


3.65(0.01) 


4.92(0.02) 


4.16(0.02) 


1.92(0.02) 


1.06(0.001) 


G345.00-0.23 


-27.69(0.04) 


5.56(0.05) 


1.53(0.02) 


1.36(0.02) 


0.46(0.02) 


0.26(0.02) 


NGC 6334 FIR-V 


-6.63(0.01) 


3.45(0.01) 


2.08(0.02) 


1.61(0.02) 


0.76(0.02) 


0.39(0.02) 


NGC 6334F 


-7.00(0.03) 


5.08(0.03) 


2.17(0.02) 


1.65(0.02) 


0.83(0.02) 


0.45(0.02) 


G35 1.77-0.54 


-3.16(0.02) 


4.88(0.02) 


3.10(0.03) 


2.56(0.02) 


1.13(0.02) 


0.62(0.02) 


G353.41-0.36 


-16.78(0.02) 


4.02(0.02) 


1.89(0.02) 


1.57(0.02) 


0.58(0.02) 


0.19(0.02) 


W28 A2(2) 


8.88(0.01) 


3.65(0.01) 


3.73(0.02) 


3.21(0.02) 


1.51(0.02) 


0.76(0.02) 


W31 (2) 


-2.78(0.03) 


5.00(0.03) 


2.08(0.02) 


1.76(0.02) 


0.81(0.02) 


0.43(0.02) 


W33 CONT 


34.80(0.01) 


3.79(0.01) 


4.20(0.02) 


3.77(0.02) 


1.77(0.02) 


1.05(0.02) 
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Table 8. Fractional SiO and CH 3 CCH abundances. The SiO column densities derived from the 29 SiO(2 - 1) line centres and the 
CH3CCH column densities derived from the J = 5k - 4 K lines are used for these estimates. The H2 column densities are derived 
from the 1.2 mm dust emission intensities assuming that T dust = T mt , where T I0t is derived from the CH 3 CCH(5^ - A K ) lines. 



Source 


Trot 


7V(H 2 ) 


iV(CH 3 CCH) 


^(CH 3 CCH) 


JV(SiO) 


*(SiO) 




[K] 


[10 22 cm- 2 ] 


[10 14 cnr 2 ] 


[IO- 9 ] 


[10 13 cnr 2 ] 


[10- 10 ] 


IRAS 12326-6245" 


34.0 


6.4 


2.0 


3.1 


0.8 


1.3 


G326.64+0.61 


29.0 


5.7 


3.4 


6.0 


2.0 


3.5 


OH328.81+0.63 


30.5 


7.4 


4.9 


6.6 


3.8 


5.1 


IRAS 15566-5304 


27.2 


4.7 


2.2 


4.7 


3.1 


6.6 


G330.95-0.19 


31.9 


8.0 


3.5 


4.4 


1.2 


1.5 


G345.01+1.8N 


34.5 


7.3 


3.5 


4.8 


2.3 


3.2 


IRAS 16562-3959 


35.9 


8.3 


15.0 


18.0 


1.6 


1.9 


G345.00-0.23 


30.6 


6.0 


4.0 


6.7 


3.3 


5.5 


NGC 6334 FIR-V 


33.6 


9.2 


5.8 


6.3 


1.8 


2.0 


NGC 6334F 


35.8 


16.3 


6.3 


3.9 


1.6 


1.0 


G35 1.77-0.54 


34.4 


13.7 


8.9 


6.5 


5.5 


4.0 


G353.41-0.36 


25.2 


13.7 


4.2 


3.1 


3.0 


2.2 


W28 A2(2) 


34.6 


7.7 


11.2 


14.6 


1.7 


2.2 


W31 (2) 


35.0 


9.2 


6.2 


6.8 


2.9 


3.2 


W33 CONT 


38.6 


16.4 


14.3 


8.7 


1.8 


1.1 



"The rotational temperature and CH 3 CCH column density derived from the J = 6 K - 5 K line are 37.2 K and 2.5 • 10 14 cm 2 , respectively. 

Table 9. Linear radii, core masses estimated from 1 .2 mm continuum (M cont ) and assuming virial equilibrium (M v ; r ) and H2 
average number densities. 



Source 


R 

[pc] 


M mnt 
[Mo] 


M viI 

[Mq] 


n(H 2 ) 
[cm- 3 ] 


IRAS 12326-6245 


0.1 


2200 


410 


1.2- 


10 7 


G326.64+0.61 


0.2 


850 


430 


6.0- 


10 5 


OH328.81+0.63 


0.2 


1200 


320 


8.4- 


10 5 


IRAS 15566-5304 


0.2 


810 


740 


5.7- 


10 5 


G330.95-0.19 


0.2 


8100 


1200 


5.7- 


10 6 


G345.01+1.8N 


0.1 


410 


550 


2.3 • 


10 6 


IRAS 16562-3959 


0.1 


400 


400 


2.2- 


10 6 


G345.00-0.23 


0.1 


960 


870 


5.4- 


10 6 


NGC 6334 FIR-V 


0.2 


500 


420 


3.5 • 


10 5 


NGC 6334F 


0.1 


830 


360 


4.7- 


10 6 


G35 1.77-0.54 


0.03 


120 


150 


2.5 ■ 


10 7 


G353.41-0.36 


0.3 


3500 


960 


7.3 • 


10 5 


W28 A2(2) 


0.1 


970 




5.5 • 


10 6 


W31 (2) 


0.7 


1600 


3500 


2.6- 


10 4 


W33 CONT 


0.3 


4300 


900 


9.0- 


10 5 



